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Site fidelity, home range behaviour and habitat utilization 
of male harlequin toads (Amphibia: Atelopus hoogmoedi ) 
from Suriname: relevant aspects for conservation breeding 

Martina Luger, Walter Hodl & Stefan Lotters 

Abstract. Recently, neotropical harlequin toads (Atelopus spp.) have undergone drastic population de¬ 
clines. Captive breeding has been suggested as a way to conserve harlequin frogs threatened by the am¬ 
phibian chytrid fungus and global warming. A major problem is that little is known with regard to the 
natural history of Atelopus species. We studied male A. hoogmoedi site fidelity, home range and habitat 
utilization during the long rainy season (April to August) at Brownsberg Nature Park, Suriname. Dur¬ 
ing this period, males and females occupy different habitats (i.e. streams versus forest). We studied male 
toads along a stream site. In total, 19 males were observed; relative abundance was 0.475 individuals per 
meter. Fourteen males were recaptured at the same site or its vicinity at least once and 12 maintained 
home ranges (38.1 ± 17.7 m 2 in size) for up to four months. However, male-male interaction was not ob¬ 
served. In the vicinity of streams, male A. hoogmoedi were perched 51.95 ± 36.7 cm above ground in dense 
understory. In conservation breeding efforts, we recommend to keep sexes separately prior to mating 
and to simulate the rainy season. Because males occupy home ranges at the breeding site, reproductive 
success in captivity may be higher when females are moved to male containers (which should have some 
running water imitating a stream habitat). Low male density and home range size should be considered 
when populating tanks with males. Tanks for males and breeding should address that wild specimens 
were mostly found on low vegetation. 

Key words. Anura, Bufonidae, ex situ conservation, habitat utilization, home range, site fidelity. 


Introduction 

Neotropical harlequin toads, genus Atelo¬ 
pus Dumeril & Bibron, 1841 (Bufonidae), 
comprise a group of more than 100 diurnal, 
slow-moving species, which aggregate along 
streams for reproduction (Lotters 1996). 
For almost two decades now, harlequin 
toads have dramatically declined through¬ 
out their geographic range (La Marca et al. 
2005). According to the IUCN/Conserva- 
tion International/NatureServe ‘Global Am¬ 
phibian Assessment’, the genus is threatened 
with extinction with 81% of the species fall¬ 
ing under the IUCN Red List category Criti¬ 
cally Endangered (Lotters 2007). Reasons 
for these declines remain unknown in most 
cases, leading Stuart et al. (2004) to use the 
term ‘enigmatic’ declines. Among others, the 
amphibian chytrid fungus ( Batrachochytri- 


um dendrobatidis) and anthropogenic cli¬ 
mate change have been suggested to play a 
role (e.g. La Marca et al. 2005, Pounds et 
al. 2006); however, others have disputed the 
link between climate change and outbreaks 
of chytridiomycosis (Lips et al. 2008). The 
few harlequin toads with apparently stable 
populations (although in IUCN Red List cat¬ 
egories Vulnerable or Endangered) include 
those from eastern Central America, the low¬ 
er eastern Andes of southern Peru and adja¬ 
cent Bolivia, the Amazon basin and the Guia- 
nas (http://www.iucnredlist.org; last accessed 
15 May 2009). 

Because some populations which had 
been stable until recently are currently un¬ 
dergoing dramatic declines, such as the Pan¬ 
amanian Atelopus limosus Ibanez, Jaramil- 
lo & Solis, 1995 and A. zeteki Dunn, 1933 
(Lips et al. 2006, K. Zippel pers. comm.), we 
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fear that soon the entire genus Atelopus will 
become extinct (La Marca et al. 2005, Lot- 
ters 2007). As advocated in the IUCN Am¬ 
phibian Conservation Action Plan (ACAP), 
one measure to prevent amphibian diversity 
from irreversible loss is conservation breed¬ 
ing (Gascon et al. 2007). According to Lot- 
ters (2007), any Atelopus species should be 
considered a candidate for ex situ conserva¬ 
tion. There is an urgent need for information 
on natural history from the wild that can be 
applied to improve captive breeding. Only a 
few studies are available (reviewed in Lot- 
ters 1996), mainly addressing spatial and 
temporal associations with a particular hab¬ 
itat. At the Monteverde cloud forest, Cos¬ 
ta Rica, Atelopus sp. (cf. varius ) aggregates 
along streams during the dry season (No- 
vember-May) and disperses more uniformly 
during the rainy season. Individuals of both 
sexes exhibit site fidelity and maintain terri¬ 
tories (Crump 1986,1988, Pounds & Crump 
1989 - each as A. varius). In A. carboneren- 
sis Rivero, 1974 from humid forests in the 
Cordillera de Merida, Venezuela, both sexes 
maintain defined home ranges throughout 
the year, except for the somewhat drier pe¬ 
riod (May to July) when individuals migrate 
to streams and reproduce (Dole & Durant 
1974 - as A. oxyrhynchus). Sexton (1958) ob¬ 
served migrations to streams in both sexes of 
A. cruciger Lichtenstein & Martens, 1856 
from coastal forests in Venezuela during the 
dry season (December to March), and pro¬ 
posed that some male A. cruciger are pro¬ 
longed residents along streams. In Panama, 
A. zeteki aggregate along streams from late 
November/early December (the beginning 
of the mating season and the transition be¬ 
tween the wet and dry seasons) to January 
(Karraker et al. 2006). Both sexes exhibit 
site fidelity (C. L. Richards pers. comm.), 
with daily movements to and from elevated 
nocturnal sleeping positions (Lindquist et 
al. 2007). Male A. pulcher Boulenger, 1882 
from Andean Peru showed site-fidelity and 
both sexes aggregated along streams dur¬ 
ing reproduction in the dry season (May to 


July, J.L. Brown pers. comm.). However, not 
all species of Atelopus follow these tempo¬ 
ral patterns, since for example A. franciscus 
Lescure, 1974 from coastal French Guiana 
(Boistel et al. 2005) and Atelopus sp. cf. spu- 
marius Cope, 1874 from near Iquitos, Peru 
(J.L. Brown pers. comm.), reproduce during 
the rainy season. These observations suggest 
that tempo-spatial behaviour varies among 
the genus Atelopus, which is relevant for con¬ 
servation breeding. 

We studied Atelopus hoogmoedi Lescure, 
1974, a species widespread on the eastern 
Guiana Shield with seemingly intact popula¬ 
tions (Rueda et al. 2005, Luger et al. 2008). 
The senior author’s preliminary observations 
revealed that male A. hoogmoedi are perma¬ 
nently found along small streams, while fe¬ 
males are encountered in forest at some dis¬ 
tance from streams. We observed females 
migrating to streams in January 2003, at the 
start of the short dry season. We studied site 
fidelity, home range size and habitat utiliza¬ 
tion in males of this species along a small 
stream in Suriname during the rainy sea¬ 
son. Further, the relevance of our findings 
for conservation breeding of A. hoogmoedi is 
discussed. 


Material and methods 

Field work 

Field work was conducted at Brownsberg 
Nature Park, Suriname (Fig. 1). Established 
in 1970, the park currently encompasses 
12,200 ha, i.e. most of Brownsberg Hill, a la- 
terite range reaching 500 m elevation. The 
area is characterized by rugged terrain with 
evergreen forested steep slopes. The climate 
is seasonally tropical and humid with annual 
mean temperature of about 26°C and annu¬ 
al bimodal rainfall of about 2,400 mm (Fitz¬ 
gerald et al. 2002). 

Atelopus hoogmoedi males are commonly 
encountered along streams in the area dur¬ 
ing visits by the first author in 2003 and 2004 
(Luger et al. 2008). We studied a population 
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Fig. 1. (A) Suriname with major river systems (thin black lines), Van Blommenstein Meer (grey), Pa¬ 
ramaribo (square) and location of Brownsberg Nature Park (circle). (B) Aspect of Brownsberg Nature 
Park with main road (bold black line) and tracks (thin black lines), headquarter (square) and creeks 
(thin grey lines with arrows); the location of the 1000 m 2 study area along the Koemboe creek where 
Atelopus hoogmoedi was studied is indicated by a circle. 


of A. hoogmoedi from 19 April to 18 August 
2004 (47 observation days) along Koemboe 
creek (Fig. 1). This small stream, bordered 
with riparian forest, is located approximately 
2 km from the main research station of the 
park (4°56’39.8” N, 55°n’47.7” W; ca. 370 m 
above sea level). The streambed is 3 m wide 
on average; water depth varies between 10 
and 80 cm. The stream gradient results in 
numerous small cascades with splash zones 
and pools below. Stream banks consist of 
sand, gravel and large boulders. Several fall¬ 
en trees are found alongside and across the 
stream. The vegetation bordering the stream 
margin was mainly composed of Heliconia 
sp. and other herbaceous plants such as ferns 
and palms. Mean air temperature during the 
study period (taken every two hours day and 
night with four Gemini data loggers) was 
23.0 ± 0.96 °C (range 20.8-27.6 °C; n = 1,500); 
total precipitation was 1,250 mm (measured 
continuously with an automatic rain gauge). 

Based the first author’s on preliminary 
observations on A. hoogmoedi at the study 
site, a 1000 m 2 plot (25 x 40 m, with the 40 
m border parallel to the stream) was estab¬ 
lished to facilitate observations. It was divid¬ 
ed into a 5 m x 5 m grid system, which cov¬ 
ered both stream banks and the forested hill 
north of the stream. Every second to fourth 


day, from 06:30-18:00 hrs, the entire plot was 
surveyed for toads (i.e. in leaf litter, on fallen 
tree trunks, rocks and vegetation up to 2 m 
above ground). Each of the 19 male individ¬ 
uals encountered was weighed with a digital 
mini scale (Kwang HWA e68-50, readability 
0.01 g) and measured (snout-urostyle length, 
SUL) to nearest 0.1 mm using a calliper when 
collected for the first time. Mean weight was 
1.4 ± 0.2 g (range 1.1-1.9 g); mean SUL was 
26.1 ± 1.3 mm (range 24.1-28.9 mm). 

For mark-recapture purposes, the dor¬ 
sal and ventral pattern of every specimen 
was photographed with a Fuji digital camera 
when encountered. Since colour patterns in 
this species are constant and differ remark¬ 
ably among individuals, individual recog¬ 
nition was possible through comparison of 
photographs. Capture locations were marked 
with labelled flags to estimate home range. 
The position within the plot was taken with a 
measuring tape relative to the grid system. In 
addition, substrate type (leaf litter, low veg¬ 
etation, rocks, logs, trees) and height above 
the ground were recorded for every capture 
event. 

To determine the relationship between 
(1) the distribution of capture locations and 
habitat structures as well as between (2) the 
distribution and size of home ranges and en- 


213 







Martina Luger et al. 


vironment, habitat structures occupied by 
male A. hoogmoedi were quantified for each 
5 m x 5 m square by counting individual 
structures, i.e. density of understory growth 
(including herbaceous plants and palms > 
50 cm in height and trees < 3 cm diameter at 
breast height), trees (> 3 cm diameter at breast 
height), fallen tree logs (> 10 cm in diameter) 
and rocks (> 50 cm in diameter). In addition, 
the slope was quantified. We chose these pa¬ 
rameters based on observations made by the 
first author at the same study during an ear¬ 
lier visit. 


Data analysis 

Site fidelity was evaluated with a recapture de¬ 
cay plot in the manner described by Brown & 
Ehrlich (1980). From the slope of the regres¬ 
sion line in the linear area (y = 1.306-0.012 x, 
whereas: x = natural log of the number of in¬ 
dividuals in residence x days), we calculated 
the mean residence rate during the study pe¬ 
riod as e m , m = -0.012 (see Warren 1992). This 
residence rate is an estimate of the percentage 
of individuals recorded on one day that stay 
in the area until the next observation. 

Home ranges were calculated as fixed ker¬ 
nels (Animal Movement 2.0 extension for 
the ESRI GIS software ArcView; Hooge & 
Eichenlaub 1997). The kernel approach is 
a contouring method that calculates differ¬ 
ent probabilities of presence for any point 
within the home range of an individual (see 
Worton 1989, Seaman & Powell 1996). 
The contour smoothing factor h was estimat¬ 
ed individually from the data set obtained in 
the field, according to the least square cross 
validation method, LSCV (e.g. Rodgers & 
Carr 1998; http://blue.lakeheadu.ca/hre). 
Home ranges were computed for eight indi¬ 
viduals (with > 10 recaptures per individual) 
with the Animal Movement 2.0 extension for 
ArcView. Outliers (maximum 5% of original 
observation points) were removed by the im¬ 
plemented function Outlier Removal accord¬ 
ing to the harmonic mean method (Dixon & 


Chapman 1980). With this method points 
can be removed by their largest harmon¬ 
ic mean value. Home range size (n = 8) was 
correlated with an index of body condition 
of the toads using a Spearman rank correla¬ 
tion. The index of body condition was cal¬ 
culated as: body weight/SUL b x 10 3 with b 
= lg[body weight]/lg [SUL] (see Hemmer & 
Kadel 1972). 

Habitat utilization: We wanted to deter¬ 
mine if correlations existed between (i) cap¬ 
ture locations and physiographic habitat pa¬ 
rameters, as well as between (ii) distribu¬ 
tion of home ranges (presence/absence) and 
physiographic habitat parameters. Therefore, 
we quantified habitat structures, including 
density and position relative to the stream 
of understory growth, trees, dead wood and 
rocks. We performed two principal compo¬ 
nents analyses on (1) habitat structure data 
from 5 m x 5 m squares occupied by Atelo- 
pus hoogmoedi and (2) on habitat data from 
all squares, whether occupied or unoccupied. 
We used a multiple linear regression analysis 
(stepwise forward) to examine which factors 
determined the distribution of occupied loca¬ 
tions, using data only from squares in which 
toads were found. Using logistic regression, 
we determined which habitat characteris¬ 
tics differentiated occupied from unoccupied 
sites. All statistical analyses were performed 
with Statistica (StatSoft). 


Results 

Site fidelity 

During our study, we in total encountered 
19 male individuals within the 1000 m 2 plot. 
Three of them were observed only once and 
two were recaptured only once. The remain¬ 
ing 14 were recaptured more than one time. 
Along the studied stream, the relative abun¬ 
dance of male toads was 0.475 per meter. Ac¬ 
cording to the recapture decay plot, male At- 
elopus hoogmoedi showed strong site fidelity 
along the stream (Fig. 2). The low decline of 
the number of recaptures in the first few days 
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Fig. 2 . Recapture decay plot for 19 Atelopus hoog¬ 
moedi males between 19 April and 18 August 
2004. 


of observation indicates little migration or 
mortality during the study period. Further¬ 
more, the high residence rate (0.988) sug¬ 
gests strong site fidelity. 

Home ranges 

Among the eight males for which home rang¬ 
es could be computed (see Fig. 3) the mean 
size of the 95% utilization probability area 
(UPA) was 38.1 ± 17.7 m 2 (range 6.4-60.7 m 2 ). 
The 95% UPA showed overlap to some extent 
among some individuals (Fig. 3), whereas the 
activity centres (50% UPA, mean size 5.4 ± 3.4 
m 2 ; range 1.2-10.6 m 2 ) almost did not overlap 
(Fig. 3). Spearman rank correlation revealed 
a significant negative correlation between 
home range size and body condition (95% 
UPA: R = -0.762, p = 0.028; 50% UPA: R = 
-0.786, p = 0.021; mean index of body condi¬ 
tion 1.03 ± 0.1). There was no correlation be¬ 
tween residence time of individuals and the 
sizes of the home ranges. 

Habitat utilization 

Home range (Fig. 3) and capture locations 
(not shown) were non-randomly distributed. 
The two factor analyses reduced the six orig¬ 
inal correlated habitat variables to two un¬ 



Fig. 3. Home ranges of eight male Atelopus hoog¬ 
moedi (numbered) calculated as fixed kernels 
within study plot (grid). Different isopleths border 
regions of equal utilization probability, i.e. the 95, 
90 and 50% volume contours. Legend: t = fallen 
tree, i = island. 

correlated factors. Concerning factor 1 both 
analyses revealed significant factor loadings 
for the variables ‘slope’, ‘distance to the stream 
and ‘density of understory growth’. With re¬ 
gard to factor 2 this was the variable ‘trees’. 
Regarding the factors which determined (i) 
the distribution of the capture locations, we 
found that they were only explained by factor 
1 (Beta = 0.434 ± 0.167, p = 0.015; Fig. 4). Ac¬ 
cording to the logistic regression model (Fig. 



(slope, distance to creek and density of understory growth) 

Fig. 4 . Multiple linear regression: correlation be¬ 
tween capture locations (between 19 April and 
18 August 2004 in a 1000 m 2 study area; see Fig. 
1) and factor 1 (Beta = 0.434 ± 0.167, p = 0.015), 
determining the distribution of capture locations; 
significant factor loadings were obtained for the 
variables slope, distance to stream and density of 
understory growth. 
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(slope, distance to creek and density of understory growth) 

Fig. 5 . Logistic regression: correlation between 
presence/absence of home ranges (95% UPA) of 
eight Atelopus hoogmoedi (between 19 April and 
18 August 2004) and factor 1 (significant factor 
loadings obtained for variables slope, distance to 
stream and density of understory growth); distri¬ 
bution of home ranges was determined by factor 
1 (chi 2 = 25.55, p < 0.001, estimated parameter 
= -3.754 ± 1.297, t(37) = -2.895, p = 0.006). One 
outlier (x = 2.19; y = 1) was removed prior to 
statistical analysis. 

5, Maximum Likelihood, n = 39 squares, one 
outlier removed before statistical analysis), 
presence of home ranges (ii) again were de¬ 
termined by factor 1 (chi 2 = 25.55, P < 0.001, 
estimated parameter = -3.754 ± 1.297, t(37) = 
-2.895, P - 0.006). Males were found between 
o and 173 cm above the ground (mean 51.95 ± 
36.7 cm, n = 226) in low vegetation (48% of 
all encounters), on trees (6%), on rocks (17%) 
and tree logs (15%). Only 14% of the observa¬ 
tions were made in leaf litter. 


Discussion 

Male Atelopus hoogmoedi showed strong site 
fidelity during our study period. Capture 
locations showed a clustered distribution. 
Males maintained activity areas in locations 
covered with dense understory vegetation 
near the stream for up to four months. We 
found a significant negative correlation be¬ 
tween home range size and body condition. 
Specimens were frequently observed at ele¬ 


vated spots up to 1.73 m above the ground in 
low vegetation (48% of the encounters). 

Site fidelity and the maintenance of a de¬ 
fined home range are considered advanta¬ 
geous since knowledge of the environment 
favours orientation and therefore make feed¬ 
ing more efficient; escape routes and shel¬ 
ter are easier to find (Wells 2007). In our 
study species, we expect that an additional 
reason may play a role. Male mate choice is 
well known in Atelopus (e.g. Lotters 1996), 
as males try to get into amplexus with con- 
specific females by clasping almost anything 
that appears suitable. Females commonly ap¬ 
pear at the breeding sites (streams) only for 
purposes of reproduction. For these reasons, 
it might be reasonable for a male to occupy 
a defined home range next to a stream well 
before females arrive. Interestingly, since our 
observations were made outside the breed¬ 
ing season (i.e. over the long rainy season), 
we conclude that mate choice may not be the 
only reason for maintaining home ranges. 

Noteworthy, it could be costly to maintain 
a large home range as evidenced by the nega¬ 
tive correlation between sizes of home ranges 
and body conditions of toads. 

Maintenance of a home range during the 
rainy season in the proximity of running wa¬ 
ter may, in some instances, even be disadvan¬ 
tageous. Sudden floods, especially in canyons 
like in the study area, can wash away Atelo¬ 
pus populations (Duellman & Trueb 1986). 
Hence, a greater perching height may be pre¬ 
ferred as a mechanism to avoid this. It should 
also allow for an increased chance to visually 
locate approaching females during the mat¬ 
ing season and to attract them by calling ac¬ 
tivity. The fact that activity centres of home 
ranges (50% UPA) almost did not overlap in¬ 
dicates that A. hoogmoedi males show territo¬ 
riality, although aggressive behaviour among 
males, as known in other Atelopus species 
(Crump 1988, Lotters 1996) was not no¬ 
ticed at all during our study. Perhaps, active 
home range defence in A. hoogmoedi plays a 
minor role (different to home range mainte¬ 
nance) outside the breeding season. 
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Our observations on site fidelity, home 
range sizes and patterns and habitat utiliza¬ 
tion in Atelopus hoogmoedi as well as the ap¬ 
parent temporal and spatial segregation of 
males (and, by implication, of females) may 
be important for both in situ (e.g. monitor¬ 
ing, population census) and ex situ (captive 
breeding) conservation programs. Ex situ 
breeding success of A. hoogmoedi at the At¬ 
lanta Botanical Garden increased after tem¬ 
poral separation of the sexes prior to mating 
(R. Gagliardo, pers. comm.). This coincides 
well with our findings. We recommend sim¬ 
ulating conditions of the rainy season when 
keeping males separated from females. Due 
to the fact that males occupy home ranges at 
the reproductive site and that this may play 
a role in mate choice, reproductive success 
in captivity may be higher when females are 
moved to males in containers that have run¬ 
ning water for egg deposition. Low male den¬ 
sity and home range size in the wild should 
be considered when populating tanks with 
males. Besides from running water, tanks 
should allow toads to climb perching sites, 
because in the wild, males were mostly found 
on low vegetation. 
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